Theoretical Physics Institute 
University of Minnesota 



On 
On 



> 
On 



TPI-MINN-93/20-T 
May 1993 



Non-perturbative production of multi-boson states 

and quantum bubbles 



A.S. Gorsky 

Institute of Theoretical and Experimental Physics, Moscow, 117259 

and 



lO ! M.B. Voloshin 

o 

Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455 

and 

P_i| Institute of Theoretical and Experimental Physics, Moscow, 117259 

Abstract 

X 

The amplitude of production of n on-mass-shell scalar bosons by a highly virtual field 
is considered in a A0 4 theory with weak coupling A and spontaneously broken symmetry. 
The amplitude of this process is known to have an n\ growth when the produced bosons 
are exactly at rest. Here it is shown that for n ^> 1/A the process goes through 'quantum 
bubbles', i.e. quantized droplets of a different vacuum phase, which are non-perturbative 
resonant states of the field <fi. The bubbles provide a form factor for the production 
amplitude, which rapidly decreases above the threshold. As a result the probability of the 
process may be heavily suppressed and may decrease with energy E as exp(— const ■ E a ), 
where the power a depends on the number of space dimensions. Also discussed are the 
quantized states of bubbles and the amplitudes of their formation and decay. 



1 Introduction 



The problem of calculating amplitudes of production of a large number n of weakly in- 
teracting bosons has received a close attention in connection with the observation!!' @ S 
that the instanton-induced amplitudes in the standard electroweak model, when calcu- 
lated to lowest orders of the perturbation theory (in the instanton background), display 
a rapid growth with energy, associated with the growing multiplicity of gauge and Higgs 
bosons in the final state. Subsequently CornwallS and Goldbergll have pointed out 
that a similar growth takes place for the amplitudes of processes, in which many bosons 
are produced by few initial particles in a simpler setting of a A0 4 theory of a scalar field 
0. The reason for this growth is that in an amplitude of a process involving large number 
n of weakly interacting bosons the smallness of the coupling constant is compensated 
by a large number of perturbation theory graphs, which typically grows as n\. This 
growth is a manifestation of the well known factorial divergence of the coefficients of 
the perturbation theoryS, and thus for n > 0(1/ X), at which n the compensation takes 
place, the perturbation theory becomes unreliable. 

In this paper we consider the amplitudes of production of n slow on-mass-shell bosons 
by a virtual field 0, '1 — > n process', in a theory of one real scalar field with spontaneously 
broken symmetry with respect to the reflection — > — 0. The Lagrangian of this theory 
in Minkowski space-time has the well known form 

c= l -(d»<pf-\^-jy , (i) 

where A is a small coupling constant and v is the vacuum expectation value of the field. 
The mass m of the bosons propagating in either of the vacua at +v or — v is m = \^2Xv. 
A number of exact results has been obtained recently related to the amplitudes of the 
1 — > n process at the threshold, i.e. when the final n particles are produced at rest. The 
sum of the tree graphs for these amplitudes was originally found explicitly0 by using 
recursion relations@ : 



(n\ 0(0) |0) = n\ (-2 V y- n (2) 

and then reproduced within a functional technique, suggested by BrownB , an extension 
of which technique will be heavily used throughout this paper. Within this extension 
the problem of calculating the threshold production amplitudes for the theory with the 
Lagrangian (0) reduces to a Euclidean space-time calculation of the quantum average 
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of the field with the kink-type boundary conditions in (Euclidean) time: — > —v 
at t — > — oo and — > +t> at £ — > +00. The tree-level expression (0) for the threshold 
production amplitudes is determined by the coefficients of expansion in powers of e mt at 
t — > — 00 of the classical kink profile 

0o(x, £) = u tanh(mt/2) , (3) 

which provides the saddle point for the classical action S[4>]. 

The first correction to the result (0) which amounts to summing one loop graphs! 



has been found and reducestLH to multiplying the lowest order expression (Q) by the 
factor (1 + c(d) n{n — 1) A), where c(d) is a coefficient depending on the number of space 
dimensions d, 

(n\(j){0)\0) tree+lloop = n\{-2v) l - n {l + c{d)n{n-l)\) . (4) 

The n 2 X behavior of the relative magnitude of the correction suggests that at large n 
the saddle point configuration in calculation of the n-th coefficient of the expansion at 
t — > —00 of the mean field, i.e. of the quantity <ft(x)e~ s , is driven far away from the 
saddle point configuration @ of the action alone. Once the correct saddle point con- 
figuration is found, the quantum fluctuations around it should produce only subleading 
in n corrections. (A similar situation in a (l+l)-dimensional toy model was discussed in 



Ref. 12 



Here we find that in the limit of large n the correct saddle point configuration in 
calculation of the threshold production amplitudes for large n is determined by dynamics 
of the surface of the inter-phase boundary ('domain wall') separating the phases with 
= —v at large negative t and with = +v at large positive t. A WKB treatment of 
this dynamics relates it to the semiclassical properties of the spherical bubbles filled with 
the phase = +v in the vacuum with = —v. At the classical level such configurations 



were studied some time agoO-3. In particular it was found that these configurations are 
reasonably long-lived: a large bubble undergoes several pulsations of its radius before 
decaying into outgoing waves. 

In what follows it will be shown that in a theory in (d + 1) dimensions (thus d 
being the number of spatial dimensions) at n ^> 1/ A the WKB result for the amplitude 
(n\ 0(0) |0) with all particles having exactly zero momenta can be interpreted as given 
by a two-stage process. The field operator 0(0) produces a bubble and then the bubble 
couples to n particles. The resulting amplitude is large: 
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(n\ (f)(0) |0) ~ n\(-2vf- n exp(b(d) n&) (5) 

with positive coefficient b(d) depending on d. This growth is in agreement with the 
general result that in this theory the threshold amplitudes should grow not slower 



than as n\. However exactly at the threshold the phase space of the final n particles 
is vanishing and to estimate the total probability of the process one needs to know 
the behavior of the amplitude above the threshold. The presence of the bubble in the 
intermediate state implies existence of a form factor, which cuts off the amplitude above 
the threshold when any of the momenta of the final particles is larger than the inverse of 
the radius of the bubble, r -1 . An estimate of the total probability of the process 1 — > n 
with the two-stage picture can be done by evaluating the probability of creation of a 
bubble. The probability of creation of a bubble with energy E evaluated by means of 
the Landau- WKB techniqueEH is found in this paper to be given by 
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| {B(E)\(j)(0)\0) | 2 ~ exp{-2b{d) (E/m)^) (6) 

with precisely the same coefficient b(d) as in eq.(^), and B(E) stands for the state of 
the bubble with energy E. 

Equation @ shows that due to the form factor provided by the bubble not only 
the growth of the probability with energy is eliminated but in fact the total probability 
associated with multi-boson states rapidly falls with energy in the non-perturbative 
asymptotic regime. Therefore we conclude that it is extremely plausible that in this 
theory the non-perturbative contribution to the processes with production of many soft 
final particles by few initial ones does not become large at high energy in spite of the 
indications to the contrary in lowest orders of the loop expansion. Whether this behavior 
is universal and is valid for other theories, in particular for the electroweak theory, is yet 
to be studied. 

The suppression of the non-perturbative processes of the type few — *> many at high 
energy does not contradict to a possible growth with energy of the probability of the 
processes many — > many, as is discussed in the concluding section. 

The rest of this paper is organized as follows. In Section 2 a brief review is given of the 
standard perturbative calculation of the amplitudes of the 1 — > n processes within the 
equivalence to the problem of calculating the quantum average of the field with the kink- 
type boundary conditions in the Euclidean space. In Section 3 the problem of finding 
the proper saddle point configuration for calculation of the amplitudes is formulated 



3 



within the so-called thin wall approximation. The search for this configuration leads 
to considering dynamics of bubbles in the Minkowski space-time. This dynamics is 
discussed within the Bohr - Sommerfield quantization in Section 4. The properties of 
bubbles are quantitatively related to the n particle production amplitudes in Section 
5 and in Section 6 the probability of creation of a bubble by a highly virtual field is 
estimated by the Landau- WKB method. Section 7 contains a discussion of the results 
and of the validity of the approximations made in present calculations. 

2 Perturbative calculation 

In this section we recapitulate the perturbative calculation of the threshold production 
amplitudes within the equivalence Gil of this problem for the theory described by the 
Lagrangian (H) to the Euclidean-space calculation of the quantum mean field with bound- 
ary conditions corresponding to a domain wall separating vacua +v and —v. Fixing for 
definiteness that the amplitudes are calculated in the 'left' vacuum, i.e. at ((f)) = —v the 
equivalent problem can be formulated as follows. 

One first calculates in the Euclidean space-time the quantum mean field 

J (j 0(x, t) d d x] e~ s W Vcj) 
&(t) = v I - (7) 

where d is the number of spatial dimensions in the problem, V is the d- dimensional spatial 
normalization volume, and the boundary condition in the path integral at t — > — oo is 
specified by the requirement that the asymptotic behavior of there is given by 

-v + ze mt + 0(e 2mt ) (8) 

with z being a constant and v and m are the renormalized v.e.v. and the boson mass. 
The production amplitudes are then givenS by the coefficients of the expansion of 
in powers of e mt at large negative t, i.e. if one writes the expansion as 

oo 

$(t) = 5> n e" m * , (9) 

ra=0 

then 

(n\<f)(0)\0)=n\c n /z n . (10) 
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The coefficient z = c% in the leading asymptotic behavior @ is thus a normalization 
of a one particle state: without dividing by the factor z n in eq.flTOD one would obtain 
the amplitudes with the particle states normalized as (1|0(O)|O) = z. As it is written 
the equation ( |T0D gives the amplitudes with the standard normalization of the particle 
states. It is also assumed here that the normalization volume V is set to unity. 

The application of the described equivalence is straightforward at the tree level, 
which corresponds to substituting for the field <j) the uniform in space (but not time) 
solution (HI) of the classical Euler-Lagrange equations for the Lagrangian (|1]). Expand me; 
this solution in powers of e mt reproduces through eq.([L0|) the tree-level result @ for the 
production amplitudes. 

The solution @ is the familiar kink profile of the inter-phase boundary ('domain 
wall') separating the vacua —v at t — > — oo and +v at t — > +00. A particular choice of 
the normalization factor z fixes the position of the domain wall in time, i.e. fixes the 
translational zero mode of the field, which corresponds to an overall shift in the time 
direction. Therefore one can either choose a particular value of z or alternatively fix 
the overall position of the field in time. Throughout this paper we choose to fix the 
position in time by requiring that the integration in the path integral in eq.(^) runs over 
field configurations such that </>(x, 0) is zero at the boundary of the normalization spatial 
bounding box. At the classical level this fixes the center of the domain wall at t = 
and thus also sets z = 2v. The quantum corrections to the mean field (0) in general 
renormalize the factor z. 

To account for the quantum effects in the mean field (0) one writes the full field <p 
as a sum of classical and quantum parts 

0(x,t) = o (x,t)+0 9 (x,t) (11) 

and evaluates the mean value of the quantum part of the field by perturbation theory in 
the background field (f>Q. At the one- loop level this calculation 0' leads to the result 
in eq.(f|). 
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3 Semiclassical configurations for quantum effects at 
large n 



As is discussed in the Introduction the growth with n of the loop corrections within 
the quantization around the saddle point configuration (H) of the action indicates that 
at large n the n-th coefficient of the expansion of $(£) in powers of e mt at t — > — oo is 
contributed by a semiclassical configuration which strongly differs from that in eq.(|^). To 
find this appropriate 'distorted' configuration we first consider the one shown in Fig.l, 
where the domain wall (the surface corresponding to 0(x, t) = 0) assumes a non-flat 
shape slowly varying with x. The surface can be described by its x-dependent deviation 
from t — 0, i.e. by the solution t = — /i(x) of the equation 0(x, t) = 0. Assume now 
that one fixes the shape of the boundary corresponding to a particular function h(x) 
and minimizes the action with respect to all other variables of the field by solving the 
Euler-Lagrange equations with the surface of zeros being fixed. Then for large negative 
t at the point x corresponding to the maximum of h(x) the n-th harmonics of the field 
is given in the leading exponential approximation by 

2v(-l) n - 1 exp (nm(t + /i(x ) + 0(nm5)) ~ 2t»(-l) n " 1 exp(nm/i(x ))e nmi . (12) 

This behavior can be understood by considering that at the point x the evolution of 
the field in time from = towards <p — ~ v proceeds over the time \t\ — /j(xq) + 5, 
where 8 reflects the uncertainty related to the curvature of the surface of the inter-phase 
boundary. This uncertainty is of a subleading importance in situations where /i(xo) is 
large, which as will be seen is the case for n ^> 1/A. 

Equation ( |P2"D tells that the coefficient c n gets multiplied by the factor exp(nm/i(xo)). 
Thus in the leading WKB approximation this coefficient can be evaluated as 

c n ~ max [exp (nm/i(x ) — S[h] + So)} , (13) 

where S[h] is the action for the field configuration described by the shape /i(x) of the 
inter-phase boundary and So = S[h = 0] is the action of the unperturbed classical 
solution ([D, and all the pre-exponential factors are omitted. 

The appearance of the action So with the plus sign in eq. fll3f) is due to the fact that 
in the expression fl7|) for the mean field the path integrals in both numerator and the 
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denominator are calculated with the kink-type boundary conditions. Thus the factor 
exp(S'o) appears from the saddle point value of the denominator. 

If h(x) varies at scale larger than the thickness of the wall, m~ x the action S[h] can 
be calculated in the thin-wall approximation^ as the surface tension of the wall /i 
times its area A: 



S[h] = fiA[h] 



(14) 



where 



-Mt)\ + 



,2\2 



dt 



2Xv i 



nT 
3A 



(15) 



: \dV^~' ) ' 4 

Finding the maximum of the expression (0) with the action ( |i~4~D is equivalent to 
solving a surface tension problem for a (^-dimensional film in [d + l)-dimensions. The 
edges of the film are fixed at the boundary of the bounding box: h(boundary) = 0, and 
at the point x the force equal to nm is applied downwards. The maximal deviation 
ho = /i(xo) of the film will be largest if the force is applied to the center of the film, 
therefore we set xo = 0. (In fact for d > 2 the equilibrium shape of the film does not 
depend on x if this point is sufficiently far from the edges. For d < 2 there is an infrared 
behavior in this problem, so that the equilibrium deviation explicitly depends on the 
size of the bounding box. Also it is explicitly assumed throughout this paper that d > 1. 
A one-dimensional 'film' lacks intrinsic curvature, which makes most of the formulas in 
this paper singular in the formal limit d — > 1, i.e. that of a (1+1) - dimensional field 
thory, for which the present analysys is thus not directly applicable.) Assuming that the 
bounding box is spherically symmetrical in the spatial d dimensions with a large radius 
R, one concludes that the shape, which the film takes under the force applied at the 
center, is also spherically symmetrical and can be characterized by the radius r(t) of its 
slice at t = const if the slice is positioned at an instant t such that — h(0) < t < 0. The 
boundary conditions for r(t) being r(— ho) = and r(0) = R. 



In terms of r(t) the quantity S[h] —nmh entering the expression (|T3|) can be written 



as 



[r] = f l d nr d ~ l Vl + r 2 dt 

J -ho 



Eh 



(16) 



where the factor nm is identified with the total energy E, r = dr/dt, and Id is the 
(d — 1) dimensional volume of unit sphere Sd-i- Id = 2 n d ^ 2 /T(d/2). The integral in 
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eq.flTBD is nothing else than the Euclidean action for a spherical bubble in the thin wall 
approximation!!!! in the theory with degenerate vacua. Since E is the conserved value 
of the Hamiltonian for the classical trajectory r(t), the functional s[r] is identified as the 
truncated action 

s[r] = J p E dr , (17) 
where Pe is the Euclidean momentum conjugate of r: 

= (18) 

On an Euclidean-space classical trajectory the value E of the Hamiltonian is related to 
Pe and r as 

^ 2 + p|=(^/ir d " 1 ) 2 . (19) 
The solution of the latter equation for pe in terms of E and r reads as 

PE = yJ(ldt*r d -i)> -E^ = l d ^^Jr^-rl d ~ 2 , (20) 

where 

r = (E/(l d im))£* . (21) 

This solution immediately reveals an important point: there is no real solution for pe 
and thus for r(t) at r < r . In the equivalent surface tension problem the origin of this 
behavior is obvious: there is a minimal radius equal to tq of a slice of the surface with 
surface tension \i that can support the force E. In terms of quantum mechanics with 
the action fllTD the point r corresponds to the classical turning point, and for r < r 
the evolution of the system proceeds in the Minkowski time. We are thus compelled 
to consider the evolution of the radius of the bubble along the complex time trajectory 
shown in Fig. 2, on which the part of the trajectory with r < r evolves along imaginary 
Euclidean, i.e. real Minkowski, time. In the Minkowski space this part of the trajectory 
describes the bubble expanding from size r = to the classical turning point r = r . 
Therefore before proceeding with further evaluating the coefficients c n by eq.fll3|) it is 
appropriate to discuss few properties of the bubbles in the Minkowski space-time. 
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4 Quantizing the bubbles 



The Hamiltonian dynamics of the bubble in the Minkowski space-time in the thin wall 
approximation is described by the simple substitution t — > it in the previous Euclidean 
space formulas. In particular the analog of the equation (|i~9"f) for the Hamiltonian is 

H 2 - P 2 M = (l d »r d - 1 f (22) 

with the Minkowski momentum pu- The classical Minkowski-space trajectory with 
energy E corresponds to oscillations of the bubble between the turning point r = tq 
(eq.(^l|)) and r = 0. Naturally it could be that instead of oscillating the bubble would 
quickly dissipate into outgoing waves. However a numerical study of the classical evo- 
lution of the field of the bubble-type configuration® (not constrained by the thin wall 
approximation) has revealed that the bubbles undergo at least several oscillations be- 
fore they emit a larger portion of their energy in outgoing waves. This implies that the 
lifetime of a bubble is at least longer than the period of oscillation T ~ r$. Therefore 
we start with discussing the bubbles as if they were stable and later take into account 
their slow decay. 

The part of the trajectory near zero radius, r < m _1 , cannot be described within 
the thin wall approximation since the thickness of the wall is of order vrC 1 . However, 
at large energy, corresponding to the turning radius tq 3> m~ l , most of the evolution of 
the bubble proceeds within the applicability of the thin wall approximation, and this is 
the part from which most of the action comes. According to equation ([H]) the condition 
tq ^> mT x implies that 

n = E/m > m 3 - d /X , (23) 

where the right hand side is the inverse of the dimensionless coupling in the theory and 
thus is assumed to be much bigger than 1. In particular for d = 3 (the normal (3+1) 
dimensional theory) the condition (|23|) translates into n ^> 1/A. Thus the applicability 
of the present calculations lies within an essentially non-perturbative domain. 

The oscillatory motion of the bubbles can be quantized and the discrete energy levels 
found by applying the Bohr - Sommerfield quantization rule: 

I(E)=<fp M dr-2Tru(E) = 2nN, (24) 
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where the integral runs over one full period of oscillation and contains the momentum pm 
determined by eq.(p2|) in the thin wall approximation. The quantity v{E) is a correction 



to the thin wall limit, which arises from the contribution to the action of the motion at 
short distances r ~ m -1 , where the latter limit is not applicable. Since at such distances 
Pm ~ E, by order of magnitude v{E) can be estimated as v{E) ~ E/m. The integral 
in eq.(p4t) is of the order of Er^, and is thus much larger than v{E) once the condition 



ro ^> m 1 is satisfied. In terms of the turning radius r$ the quantization relation (|24]) 
reads as 

k dl ir* = 2tt (N + u{E)) , (25) 
with k<i being a numerical coefficient, 



, _ j v^Fr[l/(2d-2)] 

d ~ d 2(d-l)T[3/2 + l/(2d-2)] ■ [ DJ 

At a large energy, corresponding to the condition (|23|), the correction term with v{E) 
can be neglected and one finds the expression for an energy level En in terms of the 
number N of the level: 

E N = (— ) • (27) 
By the relation (^5|) the condition ro ^> mT 1 requires that 

N > m 3 - d /X , (28) 



where the right hand side by itself is assumed to be a large number. The energy E 



N 



of the level then satisfies the condition (23). It can be readily noticed that for such N 



the spacing between consecutive levels is small in comparison with the mass m of the 
bosons: AEn ~ r^ 1 m. This, perhaps, in addition to the reflectionless property of 
the wall explains the relative stability of large bubbles. Indeed, emission of individual 
quanta would require transitions between states with a large difference of their quantum 
numbers iV: AiV ~ mr ^> 1. The overlap integral for the wave functions of the 
levels with large difference of their numbers is exponentially small in AiV. Therefore the 
probability of emission of the bosons by large bubbles should be strongly suppressed by 
the parameter mr§. However there perhaps can be a larger probability of simultaneous 
emission of many bosons at short distances, i.e. when the bubble contracts to r < mT 1 
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during the oscillations. Admittedly at present there is little understanding of the decay 
of the bubbles besides the studyO by means of classical field equations. 



5 Quantized bubbles and multi-boson processes 

Returning now to evaluation of the coefficients c n by eq.fllBD, we notice that the exponent 
there receives a real contribution from the action on the Euclidean part of the trajectory, 
i.e. when r changes between r = r$ and r = R, and an imaginary part on the Minkowski 
part of the trajectory, i.e. between r = and r = r . Let us first evaluate the Euclidean 
part, where the truncated action (eq.(|l7])) is given by 

s[ r ]= f R p E dr= f R ^{l d nr d - v ) 2 - E 2 dr . (29) 

Jro J ro 

The latter integral diverges in the limit R — > oo. However the quantity of interest is 
the difference So — s[r], which enters the exponent in eq. (|13]) . Since the action So of 
the unperturbed configuration formally corresponds to the truncated action s[r] at zero 
energy E = and has the same leading divergence, one can first calculate the derivative 
of s[r] with respect to E and then find the difference of the action as 

So ~ s[r] = [ E dE . (30) 



/o V dE r 

For the derivative with respect to energy one finds from eq 



ds[r] E f R dr 

(31) 



dE Id fj, Jr ^J r 2d-2 _ r 2d - 2 

Before discussing this integral for large energy, it can be noted that quantitatively 
this expression can be used to estimate the quantum effects in the coefficients c n also at 
low energies where formally ro is less than thickness of the wall. Then the integral should 
be cut off at a lower limit r = r\ ~ m^ 1 , which does not depend on E if E/m m 3 ~ d / A. 
Then one finds 

rm / E 2 fR dr \ 

where is the tree- level expression for the coefficient c n . When the exponent is ex- 
panded in powers of its argument this gives the (1 + const ■ n 2 \) behavior of the first 
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quantum correction to the production amplitudes, which one finds by a calculation of 
one-loop graphs (eq.(^|)). Moreover for d < 2 the integral in equations fl3~T|) and 

(|32|) still diverges in the limit R —>■ 00. Therefore the dependence on the infrared cut 
off R dominates the quantum effects and makes them insensitive to the region of small 
r. In this case in the leading at large R approximation the exponent in eq.flTBD can be 
found for any E with the obvious result 

c " ^ = 4 exp [ e (2 - 1) u m ) ■ (33) 

It is also clear that the Minkowski part of the trajectory r(i) contains no infrared de- 
pendence and thus does not contribute to this result in the leading in R approximation. 
It can be readily verified by the technique of Refs. |TT)| , [TT| ] that the infrared behavior of 
the one-loop corrections to the coefficients c n exactly reproduces the expansion in the 
latter equation of the exponent up the first power of its argument. 

Therefore considering the case where the equation (^Tj) is applicable down to r with 
the condition that r ^> m~ x is sensible only when the integral in that equation is finite 
in the limit R — > 00, i.e. when d > 2, which includes the most interesting case of d = 3. 
Setting R = 00 in eq.([H]) one finds for d > 2 

with the dimensionless factor fd given by 

0Fr[l/2-l/(2d-2)] 

Jd 2(d-i)r[i-i/(2d-2)] ■ { } 

Therefore one finds the real part of the difference Sq — s[r] associated with the 
Euclidean part of the trajectory r(t) to be given by 

&--W' = '< i -T LE {£r ■ < 36 > 

Let us now evaluate the contribution to the coefficients c n of the Minkowski part of 
the evolution of the bubble. The simplest trajectory, which links r = with r = ro 
consists of one half of the period which contributes to the coefficient c n the factor 

F = exp(z/(£)/2) , (37) 
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where the action integral 1(E) over a full period is defined in eq.fl24|) However there are 
also trajectories connecting r = with r = r , which differ from this one by an integer 
number k of full periods of oscillation. Each of these trajectories provides a saddle point. 
Therefore one has to take the sum of contributions of these saddle points. The sum has 
the form 

^■f^lHlfflW- ^ (38) 

It is clear that this expression sums to infinity when the Bohr-Sommerfield relation 
(eq.(§4])) is satisfied, i.e. 1(E) = 2nN. In other words the factor F(E) develops poles 
at the values of energy coinciding with the positions of the bubble levels En given by 



eq. ( 27 ) . It can be reminded that thus far the decay of the bubble levels is completely 
ignored. Therefore we come to the conclusion that in this approximation the amplitude 
of production of n static bosons consists of poles at the energies E = En with the 
residues proportional to (— 1)^ exp(S , — s[r])#, the latter being given by the equation 
flSHp. The sign alternating factor (— 1)^ arises from exp(iI(Eff)/2) = exp(inN) in the 
numerator of eq. (|38|) . The decay width of the bubble levels can be taken into account in 
the Breit-Wigner approximation by shifting the positions of the poles into the complex 
plane: En — > En — iTn /2. 

To conclude this section we collect all the factors in our estimate of the production 
amplitudes of n bosons, all being at rest, at energy E corresponding to n = i?/m > 
m 3 ~ d / A and write the final result in the form 



^ (1 E (yld/Ji) 



(39) 



This expression consists of poles corresponding to the energy levels En of quantized 
bubble, i.e. corresponding to 1(E) = 2ttN. Therefore one can conclude that the pro- 
duction of the bosons in ultra-high-energy limit goes through intermediate states, which 
are the quantum levels of a bubble. 

6 Amplitudes of bubble formation and decay 

The residue of an individual pole in eq. (|39|) is given by 
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res N = (-If exp(D(E N )) , (40) 

where the shorthand notation D(E) is used for the expression in eq. (|36|) , which also 
enters the exponent in eq.([39|). The residue can in general be written as 

resiv = A(l -> B N ) ■ A(B N -> n) , (41) 

where A(l — > £>at) = (-8at|0(O)|O) is the amplitude of production of the Af-th state of 
the bubble, Bn by a virtual field and A(Bn — > n) is the amplitude of transition of this 
state into n bosons, all being at rest. The former amplitude can be evaluated by the 
Landau- WKB technique^' for calculating transition matrix elements. According to this 
technique in the exponential approximation the matrix element of an operator / between 
states with energy E\ = and E 2 = E, (i£|/|0), is found by matching the Euclidean 
classical trajectories, one with the energy of the initial state, i.e. E\ = 0, and the other 
with the energy of the final state E 2 = E, which runs between the matching point and 
the turning point. The matrix element is given in the exponential approximation by 

(£|/|0)~exp(s(£)- S (0)) , (42) 

where s(E) is the Euclidean-space truncated action on the trajectory. The specific form 
of the operator / enters only the pre-exponential factor (unless the operator / itself is 
exponential, which is not the case in the problem under discussion) and can be ignored. 

The configuration shown in Fig. 3 displays such matching of the evolution in the Eu- 
clidean space. The evolution starts with zero energy, which corresponds to a flat domain 
wall. The field then matches that of a bubble with energy E, whose radius then con- 
tracts down to the turning point tq. The configuration is then symmetrically extended 
beyond the turning point, so that it gives the square of the matrix element. Clearly, 
the difference of the truncated action (s(E) — s(0)) on one half of this configuration, i.e. 
from t = — oo to the turning point, is equal to minus that given by equation d36|) . Thus 
one immediately finds the estimate 

\A(1 - B N )\ = 1(^10(0)10)1 ~ exp(-D(E N )) . (43) 

A comparison of the formulas (|40"D, ( fl~l]) and (fHf) leads to the following reasoning. 
The product ( fill) of the amplitudes of formation and decay of the bubble is exponentially 
large in E d ^ d ~ 1 ^ (eq.(|40|)) while the formation amplitude given by eq. (|43|) is exponentially 
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small. Thus the amplitude of the coupling of the bubble to n bosons, which are all being 
at rest contains the doubled positive exponent: 



To reconcile this extremely strong coupling of the bubble to the state of n bosons, 
in which they all have exactly zero spatial momenta, with a total decay rate that is 
not exponentially large, one inevitably has to assume that the coupling of the bubble 
to bosons develops a form factor which sharply decreases above the threshold and thus 
is capable of suppressing the double-exponential and the factorial energy growth of the 
amplitude at the threshold. In view of this observation it is extremely likely that in the 
processes 1 — ► n, whose amplitude at the threshold has only single exponential growth 
factor, the same form factor makes the total probability exponentially suppressed at 
high energy. 

7 Discussion and conclusions 

The results of the search for the correct saddle point for calculation of the n-th coefficient 
c n in the expansion of the mean field <fi under the kink type boundary conditions justify 
the thin wall approximation used in this paper in the limits considered. Indeed, the 
uncertainty of this approximation, expressed by 5 in eq . (|TTD , which is of the order of the 
thickness of the wall, becomes small in comparison with the main term if the maximal 
deviation of the wall ho is much larger than m -1 . One can readily see that this is indeed 
the case in the calculations of this paper. For the number of space dimensions d < 2 the 
maximal deviation grows with the size R of the bounding box and thus for the leading 
infrared terms the thin wall approximation is applicable at any n. As mentioned in 
connection with the equation ( |33"D the result for this case can be checked against a direct 
calculation©' of the infrared behavior of the amplitudes at the one-loop level. In 
the infrared finite case of d > 2 in particular for d = 3 the applicability of the thin wall 
approximation is guaranteed in the limit of large n: n ^> m 3 ~ d /\. This follows from 
that the maximal deviation ho is of the order of vq, thus according to equation fl2l|) 



There is possibly not a coincident correlation between the qualitative dependence of 
the present calculations on the number of spatial dimensions and a similar dependence 
of the spectra of states of many soft bosons with point-like attraction. Noticing in this 




(44) 



ho m ~ (n\/m 



3-dU/(d-l) 
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connection that the interaction between soft bosons near the threshold in the theory with 
the Lagrangian ([I]) is attractive, we recall that for d < 2 an attraction at short distances 
produces bound states for any number of bosons, starting from n — 2. The existence 
of such bound states strongly changes the spectral density of the n-particle scattering 
states, which is another way to formulate that there is a strong behavior of the form factor 
for production amplitudes near the threshold. For d > 2, in particular for d — 3, a weak 
attraction at short distances is insufficient to bind few bosons. However if the number of 
bosons is non-perturbatively large, n ^> m 3 ~ d /\ they can form collective multi-particle 
bound states, which in the terminology of the present paper are the bubbles of the field 
<fi. Though a detailed relation to the present calculations is yet to be understood, this 
may provide another explanation of why the semiclassical analysis can be applied for 
any n if d < 2 and why it becomes applicable in the case of d > 2 only when the number 
of produced bosons is non-perturbatively large, n ^> 1/A. 

The main result of the present paper is the following hierarchy of the amplitudes for 
transitions between one highly virtual particle, the state of n bosons, all having zero 
spatial momenta, and a state of a spherical bubble with energy E = nm: 

A(l -> n) ~ A(l -> B) ■ A(B -> n) ~ e D{E) , 

A(B -> n) ~ e 2D{E) . (45) 

When the particles have non-zero momenta a form factor arises, due to the size tq of the 
bubble which cuts off the phase space integration. The suppression due to the form factor 
can be evaluated from the fact that the decay rate of a large bubble is not exponentially 
large in its energy, i.e. in this exponential scale the rate is T(B — > many) ~ 0(1). If 
the processes 1 — > many at ultra-high energy are going through the bubbles, as it is 
strongly indicated by the calculations in this paper, their rate should be cut off by the 
same form factor, which implies 

T(l -> many) ~ \A(1 -> B)\ 2 T(B -> many) ~ e ~ 2D{E) (46) 

and thus these processes are strongly suppressed. 

This hierarchy can be extended to processes many —>■ many, which potentially can 
go at high temperature. If such scattering processes are also mediated by the bubbles, 
one can estimate 
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A(n -f n) ~ A(n -> 5) • A(S -> n) ~ e 4D(jB) . (47) 

A kinetical calculation of the rate in a thermal equilibrium in this case involves the form 
factor for both the final and the initial states. Thus the rate of the multi-particle pro- 
cesses at high temperature should be T(many — *> many) ~ 0(1), which is in agreement 
with standard thermodynamical calculations. 

For a very approximate understanding of the role of the form factor one can invoke 
the following reasoning. The amplitude of the process 1 — > n above the threshold is 
a function of all the n spatial momenta of the final bosons, A(pi), and the total 
probability is given by 

1 r n H d n- 

r = i [ /|A(p,)| 2 n 7 4f^ ( 48 ) 

n\ J 2ei (2vr) d 

with the amplitude A(0) corresponding to all the momenta vanishing being given by 
eq. fl39|) . The presence of the bubble with the radius Tq in the intermediate state makes 
the amplitude A(p») to sharply decrease when any of the momenta is of the order of 
1/ro- Therefore the phase space corresponding to this region in the momentum of each 
particle behaves parametrically as (r )~ nd ~ E^ Ed ^ d ~ 1 \ This suppression is sufficient 
to eliminate the factorial growth of the amplitude A(0), but is not sufficient to overcome 
the exp(— const ■ E d ^ d ~~ 1 * > ) enhancement. The estimate of the total rate in eq.(fSJ) shows 
that it is most likely that the actual suppression due to the finite size of the bubble is 
somewhat stronger than in this simplistic reasoning, perhaps, due to coherence effects. 

As is mentioned before our estimate of the form factor suppression relies on that 
the decay rate of a bubble does not grow exponentially as some power of the energy 
of the bubble. Therefore it is appropriate to emphasize once again the arguments in 



favor of this behavior. One argument is based on the numerical observation!^ of the 
classical relative stability of the bubble: it dissipates its energy over several oscillations, 
thus its lifetime is at least not shorter than of the order of the period of oscillations 
T ~ r . In the quantum theory the suppression of emission of the bosons from large 
distances follows from the fact, discussed in Section 4, that the splitting of the levels of 
the bubble at high energy is much smaller than the mass m of the quantum. Therefore 
emission of a boson requires a transition between levels with a large difference AiV of 
their quantum numbers, which is exponentially suppressed in AiV. Thus the emission 
can take place during the time when the radius of the bubble is small r < 0(m~ l ) and 
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the walls collide into each other. However even if the probability of annihilation of the 
walls into outgoing bosons in this region is of order one, the bubble spends only a small 
fraction of time at such small radius, so that the total decay rate is proportional to the 
inverse of the period T. Therefore there seems to be every reason to assume that the 
decay rate of a bubble is limited by const /T. 
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Figure captions 

Figure 1. A configuration of the field, corresponding to the inter-phase boundary bent 
to the maximal deviation h . The evolution of the field between the point at a negative 
time t (heavy dot) and the boundary proceeds over the time \t\ — h + 5, where 5 is of 
the order of the thickness of the domain wall. The bold vertical lines at the edges are 
the world lines of the boundaries of the spatial bounding box. 

Figure 2. The world line of the bubble walls in a complexified space-time. The bubble 
evolves in the Euclidean space-time at r > r and in the Minkowski one, when r < r . 

Figure 3. The Euclidean field configuration for calculating the square of the matrix 
element (B\(f)\0) by the Landau- WKB formula. At large negative time the field evolves 
by the classical solution with zero energy, corresponding to a flat domain wall (lower ho- 
risontal line). Then it matches on the configuration with a large energy E, corresponding 
to a bubble, which contracts down to the turning radius tq. Beyond the turning point 
the configuration is symmetrically reflected in time, hence it represents the square of 
the matrix element. The plus and minus signs indicate the phases, corresponding to the 
field approaching +v or —v. 
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